Dynamics of magnetization in frustrated spin-chain systems: Ca 3 Co 2 06 
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The magnetization dynamics of the triangular lattice of Ising spin chains is investigated in the 
framework of a two-dimensional model. The rigid chains are assumed to interact with the nearest 
neighboring chains, an external magnetic field, and a heat reservoir that causes the chains to change 
their states randomly with time. A probability of a single spin-flip process is assumed in a Glauber- 
like form. This technique allows describing properly the steps in the magnetization curves observed 
in Ca3Co2 06 and their dependence on a magnetic field sweep rate and temperature. A transition 
from a low-temperature to high-temperature phase is also observed. 

PACS numbers: 75.25.+Z, 75.30.Kz, 75.50.Ee 



Ising spin-chains packed into a two-dimensional (2D) 
frustrated lattice demonstrate a complex magnetic be- 
havior due to a combination of a low dimensionality 
and frustration. There are few groups of compounds in 
which the triangular lattice is formed by antiferromag- 
netic (AFM) Ising spin chains, e.g. CsCoCl 3 , CsCoBr 3 
[H, or ferromagnetic (FM) ones like Ca 3 Co 2 6 HQ. A 
spin-chain system SrsRhjOis discovered recently has a 
complex magnetic structure of chains Q • An interaction 
between the nearest neighboring chains in the lattice is of 
the AFM type and much weaker than the intrachain one. 
However, it causes the frustration and a large variety of 
magnetic structures p], [H, H, 0, H, @] . 

A step-like magnetization curve in CEI3C02O6 has 
drawn recently considerable attention d, HI, H, 0, H, H, 

IE 

[111 [H, [Hf . The number of the steps depends on a 
swee p ra te of the external magnetic field and temperature 
(3.IRll0j|. Two steps become apparent in the temperature 
range from 12 K to 24 K [Io| . The first step takes place 
at the zero magnetic field. Then a plateau at about 1/3 
of the full magnetization stretches up to the magnetic 
field of 3.6 T where the second step to the saturated FM 
state occurs. At least four equidistant steps are clearly 
visible below 12 K at a moderate magnetic field sweep 
rate [H, [|| . At an extremely low sweep rate the magneti- 
zation curve becomes close to the two-step shape, similar 
to that is observed at the high temperatures [8|. A re- 
sponse to alternating magnetic fields and its dependence 
on the temperature and magnetic field were investigated 
carefully in Ref.Q- The experimental results show that 
there exist two characteristic time scales of the magne- 
tization dynamics. The first is of order of 1 s and the 
second reaches a few hours. 

A crystal structure of Ca3Co206 consists of C02O6 
chains running along the c axis. The Ca ions are sit- 
uated between them and are not involved in magnetic 
interactions. The chains are made up of alternating face- 
sharing CoC>6 trigonal prisms and CoC>6 octahedra. The 



crystalline electric field splits the energy level of Co 3+ 
ions into the high-spin [S = 2) and low-spin {S = 0) 
states. The chains form triangular lattice in the ab plane 
that is perpendicular to the chains. An in-chain exchange 
interaction between high-spin cobalt ions through the oc- 
tahedra with low-spin cobalt ions is ferromagnetic. The 
parameter of the FM in-chain coupling (Jf) was found 
from the magnetic susceptibility at high temperatures 
[IH, specific heat 0], and theoretical calculations [T^. 
These estimations are in a reasonable agreement with 
each other (Jp rj 25K). The weak AFM interchain inter- 
action causes the partially disordered AFM (PDAFM) 
structure of chains (or honeycomb magnetic structure) 
below X'ci=24 K. At T^2 ~ 12 K another magnetic tran- 
sition takes place 0, Ujj. It should be mentioned that 
the topology of the magnetic net in Ca 3 Co20fi is rather 
complex, e.g. there exist helical paths [12| . 

Taking into account the strong intrachain interaction 
one can assume the chains to be in two ordered states 
(spin- up or spin-down) at low temperatures. This is a 
basic assumption of a rigid-chain model Then the 
problem is reduced to the AFM triangular Ising model 



(1) 



where <7j = ±1 is the c-axis projection of the z-th chain 
spin, J > is the parameter of the AFM interchain cou- 
pling, B is the external magnetic field, < ij > denotes 
the summation over all the nearest-neighbor pairs on the 
triangular lattice. 

A second assumption of the rigid-chain model is that 
even at a very low magnetic field sweep rate the system is 
out of equilibrium, i.e. it is in a metastable state rather 
than in the ground state [5(. The large hysteresis loop 
and strong dependence of magnetization curve on the 
magnetic field sweep rate favor this point. Conditions of 
the metastability of the system can be formulated in the 
following form [5[ 
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Uihi < 



(2) 



FIG. 1: Fragments of the triangular lattice during the simulation at T — 4 K and 1 T/min (from left to right): the initial state 
B = T, B = 0.768 T, B = 1.824 T, and B = 2.784 T. The white and black circles stand for spin-up and spin-down chains. 



where hi = Jj2<ij> a j — B is the effective field for the 
i-th chain. This inequality should be satisfied for all 
the chains. Following the single- flip technique for the 
nonequilibrium Ising model [la] the probability Wi of a 
spin-flip event at the i-th chain was formulated as 



Wt = 



if a l h l < 0, 

1 if (Tihi > 0. 



(3) 



The rigid-chain model describes properly the four-step 
magnetization curve in CaaCc^Og. On the other hand 
one should keep in mind that once the system occurs 
in a metastable state it will dwell in this state for an 
arbitrarily long time, i.e. this is a frozen metastability 
approximation. 

Recently, static magnetization curves in Ca3Co206 
were investigated by means of Monte Carlo technique 
[l6| . It was observed that the perfect triangular lattice of 
the rigid spins produces the two-step magnetization curve 
even at low temperatures. To overcome this problem the 
nearest-neighbor interactions were radomized. The mag- 
netization curves demonstrated the four-step behavior on 
such imperfect lattice. However, it should be pointed out 
that these two types of the magnetization curves were 
observed experimentally on the same sample at different 
sweep rates Q. Therefore the four-step magnetization 
curve can not originate from the lattice imperfection. 

In the present paper, we generalize the rigid-chain 
model to include dynamical aspects of the problem. Fol- 
lowing to the Glauber theory [13] we assume that the 
chains interact not only with the nearest neighbors and 
external magnetic field but also with a heat reservoir. In 
this case the probability of a spin flip of the i-th chain 
per time unit can be written down as 



Wi 



tanh 



J 

kf 



a i 




(4) 



where a is the constant of the interaction of a chain with 
the heat reservoir, k is the Boltzmann constant, T is the 
temperature, \i is the chain magnetization. It should 
be mentioned that the ratio J/fi do not depend on the 
chain length. It is determined by the plateau length in 
the magnetization curve (AB — 1.2 T for CaaC^Og). 



In effect a temperature behavior of the system is related 
to the effective temperature kT/N where N is the length 
of the chain (number of magnetic ions). 

One can see that if T — -> the expression (0| goes 
to the function ([3]), i.e. the model of Ref.[5[ is the low- 
temperature limit of the present one. In the opposite case 
T — > oo the probability of a spin-flip event is independent 
on the magnetic field and interchain coupling. At finite 
temperatures the probability depends on the ratio of the 
effective field to the temperature and is nonzero even if 
(Xihi < 0. That is why, an evolution of metastable states 
exists within this model. 

We have performed a numerical simulation of the mag- 
netization evolution on a rhombic 96 x 96 supercell of 
the triangular lattice with periodic boundary conditions. 
Few periodic structures can appear during the simula- 
tion. To avoid artificial domain boundaries the length 
of the supercell should be multiple of a periodicity of 
these structures. In order to satisfy this condition for any 
three- or two-sublattice ordering the length was chosen 
to be multiple of 6. Preliminary calculations on 24 x 24, 
48 x 48 , and 192 x 192 supercells showed a good conver- 
gence of the results. The best coincidence of the theoret- 
ical and experimental data occurred at a = 0.21 s _1 and 
N = 30 which are reasonable values. They were used in 
all further calculations. 

The simulation started with preparation of an initial 
state. As it was discussed in Ref.|5| the PDAFM struc- 
ture is rather close to the ground state of the AFM trian- 
gular Ising model. That is why, the initial state was ob- 
tained by an evolution of the PDAFM phase at the zero 
magnetic field. When the duration of this preliminary 
evolution was more than 5a _1 , results of the simulation 
became practically independent on the initial state. The 
typical duration of the preliminary evolution in our cal- 
culations was about 100a _1 . A fragment of the initial 
state is shown in FigfT^,. It should be mentioned that 
one can prepare the initial state from an arbitrary ran- 
dom state. However it requires a longer period to get the 
initial state. 

After obtaining the initial state the magnetic field was 
switched on at a constant sweep rate. We adjusted the 
magnetic field sweep rates to the experimental values of 
Ref. Q. A typical number of time steps was about 2 • 
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FIG. 2: The magnetization curves in the increasing magnetic field at different sweep rates and T — 4K (left) and at different 
temperatures with the constant magnetic field field sweep rate of 0.1 T/min (right). 



10 6 per a full cycle (increasing-field and decreasing-field 
branches) . 

Fig. [2] demonstrates an average magnetic moment as a 
function of magnetic field M{B) calculated for different 
magnetic field sweep rates at a constant temperature (in 
the left panel) and for different temperatures at a con- 
stant sweep rate (in the right panel). One can see that 
the second and third steps in the magnetization curve 
vanish with decreasing of the magnetic field sweep rate 
or with increasing the temperature. This general ten- 
dency and shape of the magnetization curves are in a 
very good agreement with the experimental observations 

To clarify an origin of such behavior we present spec- 
imens of configurations at different magnetic fields cor- 
responding to the plateaus in the magnetization curves 
at the highest magnetic field sweep rate (see FigllJ. One 
can clearly see domains in the figures. The ferrimagnetic 
state is formed by three sublattice. It is obviously three- 
fold degenerate. While the magnetic field sweep rate is 
sufficiently fast the ferrimagnetic phase starts growing at 
a number of nucleation centers. This causes the domain 
formation. With a decrease of the sweep rate the do- 
mains expand and for the extremely slow sweep rate the 
domain walls almost disappear. It is easy to see that the 
single-domain ferrimagnetic state leads to the two-step 
magnetization curve because all the spin-down chains are 
encircled by six spin- up chains. The dotted curve in the 
left panel of FigJS] corresponding to the extremely low 
magnetic field sweep rate is similar to that was obtained 
by the Monte Carlo technique [l6l |. 

Domain boundaries contain chains with other environ- 
ments and cause appearance of the two additional steps 
in the magnetization curve. For instance, at the first 
plateau (FigJTJa) there are spin-down chains with 2 spin- 
down and 4 spin-up nearest neighbors. Their critical 
spin- flip field is 1.2 T. At the second plateau (FigOJj) 



there exist spin-down chains with 1 spin-down and 5 
spin-up nearest neighbors. Their critical spin- flip field 
is 2.4 T. The higher the magnetic field sweep rate is, the 
smaller the domains become and the more apparent the 
additional magnetization steps are. 

Here we should mention that the chains in the domain 
boundaries as well as the internal chains are in metastable 
states at the plateaus, that is, they are oriented along the 
effective fields. The spin-flip probability in the form of 
Eq.((4]) allows some fraction of exited states that, in their 
turn, leads to a creep of the domain wall. This is a new 
important ingredient which appears in the present theory. 
It explains the existence of two characteristic times in 
Ca3Co2 06- The first one is related to a single spin- flip 
evens and specifies the fast processes of about 1 s. The 
second is due to the creep of the domain boundaries (up 
to 10 4 s). 

As it was argued in Ref . [|| the transition from the low- 
temperature to high-temperature regime stems from dis- 
ordering of a fraction of the chains. This mechanism is 
essentially three-dimensional. That is why, the obser- 
vation of this transition within the present 2D model is 
surprising. The present interpretation of the transition 
differs from that was assumed in Ref.[(|. It arises due 
to a drastic increase of a domain boundary mobility at 
high temperatures that leads to establishing of the single- 
domain ferrimagnetic state during a short time. 

Another interesting features that were observed ex- 
perimentally are crossings of field-increasing and field- 
decreasing branches of the magnetization curve [![. As 
one can see in Fig |3] they are also reproduced by the sim- 
ulations. 

In conclusion, in the present paper we go beyond the 
frozen metastability approximation of Ref. [5( . The simu- 
lation of the nonequilibrium evolution was performed by 
means of a Glauber-like form of the spin-flip probability 
that allows us investigating the dependence of the mag- 
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netization curves on the temperature and magnetic field 
sweep rate in a good agreement with the experimental 
data. At the extremely low magnetic field sweep rate the 
results of the 2D simulation are similar to that was ob- 
tained by the Monte Carlo technique [l6j]. The present 
model gives reasonable values of both the critical tem- 
peratures: T C2 w 10 K and T C \ w 20 K. Although it 
is not valid at high temperatures (T > 15 K) since the 
intrachain disorder becomes essential jy] . An experimen- 
tal feature that escapes the 2D simulation is a drastic 



increase of the hysteresis loop at a very low tempera- 
ture (T = 2 K) as compared to moderate temperatures 
(T ~ 5 K). Most probably it is related to an intrachain 
dynamics of the magnetization and requires a 3D simu- 
lation to be describe properly. 
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